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The Mars Entry Atmospheric Data System is a part of the Mars Science Laboratory, 
Entry, Descent, and Landing Instrumentation project. These sensors are a system of seven 
pressure transducers linked to ports on the entry vehicle forebody to record the pres- 
sure distribution during atmospheric entry. These measured surface pressures are used to 
generate estimates of atmospheric quantities based on modeled surface pressure distribu- 
tions. Specifically, angle of attack, angle of sideslip, dynamic pressure, Mach number, and 
freestream atmospheric properties are reconstructed from the measured pressures. Such 
data allows for the aerodynamics to become decoupled from the assumed atmospheric prop- 
erties, allowing for enhanced trajectory reconstruction and performance analysis as well 
as an aerodynamic reconstruction, which has not been possible in past Mars entry recon- 
structions. This paper provides details of the data processing algorithms that are utilized 
for this purpose. The data processing algorithms include two approaches that have com- 
monly been utilized in past planetary entry trajectory reconstruction, and a new approach 
for this application that makes use of the pressure measurements. The paper describes 
assessments of data quality and preprocessing, and results of the flight data reduction from 
atmospheric entry, which occurred on August 5th, 2012. 

I. Introduction 

On August 5th 2012, the Mars Science Laboratory (MSL) 1 entry vehicle (EV) successfully entered the 
Mars atmosphere and landed the Curiosity rover safely on the surface of the planet in Gale crater. MSL 
carried with it a unique instrumentation package designed to measure the aerodynamic and aerothermal 
environments during atmospheric entry. This instrumentation package is known as the MSL Entry, Descent, 
and Landing Instrumentation (MEDLI), 2 which consists of two major subsystems, the Mars Entry Atmo- 
spheric Data System (MEADS) and the MEDLI Integrated Sensor Plugs (MISP). The MEADS consists of 
seven pressure transducers connected to flush orifices in the heat shield forebody to measure the pressure 
distribution. The MISP devices are a system of thermocouple and recession sensors to provide aerothermal 
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measurements of the heat shield performance. The MEDLI sensors provide key measurements that can 
be used for trajectory reconstruction and engineering validation of aerodynamic, atmospheric, and thermal 
protection system models in addition to Earth-based systems testing procedures. Such validation directly 
benefits future planetary entry, descent, and landing (EDL) by reducing uncertainties associated with these 
models and procedures. 

The basic MEADS science objectives are to reconstruct atmospheric data variables from pressure mea- 
surements alone when the free stream dynamic pressure is above 850 Pa. In particular, the objectives are 
to estimate the angles of attack and sideslip to within 0.5 deg and the dynamic pressure to within 2%, in 
a 3cr sense. Secondary objectives are to estimate the Mach number, freestream density and atmospheric 
winds from the MEADS pressure measurements, when combined with the on-board Inertial Measurement 
Unit (IMU) data. These measurements serve to enhance the MSL trajectory reconstruction and performance 
analysis, and enable a separation of the aerodynamics from the atmosphere, which prior to MEADS has not 
been achievable for Mars EDL reconstruction. The concept is an implementation of the Flush Air Data 
System (FADS). 

The concept of the Flush Air Data System (FADS) was conceived and developed specifically to provide 
research quality air data during the hypersonic flight regime where the classical Pitot static probe could not 
survive. While limited hypersonic air data was acquired during the X-15 program using a servoed Q-Ball 
Air Data System (see References 3 and 4), such a concept is not compatible with blunt entry configurations. 
A blunt-body FADS was attempted on the Viking 1 and 2 entry vehicles, with limited success. 5 The 
initial flight FADS was proposed and developed under the Orbiter Experiments Program as the Shuttle 
Entry Air Data System (SEADS), documented in References 6-22. The Space Shuttle Orbiter provides a 
model for the accomplishment of a staged program for verifying EDL technology readiness. In addition 
to an elaborate ground test program and sub-scale flight testing, the initial manned full scale testing was 
accomplished using full scale un-powered final approach and landing tests. To support the reentry flight, an 
expanded data acquisition program was defined for the initial orbital flights of the Orbiter. These flights 
were flown on relatively low risk entry trajectories while EDL data was collected by the Development Flight 
Instrumentation (DFI). These first five development flights were specifically designed to collect the EDL data 
required to expand the flight envelope as uncertainties were addressed and reduced. The data was oriented 
towards the validation of structural loads, heat loads and GN&C uncertainties. The data from these flights 
significantly impacted the orbiter EDL operations. As the design and development of SEADS hardware and 
data reduction software evolved the system’s objective expanded to demonstrate the concept as one with an 
across the speed range air data capability - subsonic through hypersonic and ascent. SEADS successfully flew 
five research/demonstration of concept missions on the orbiter OV-102 (Columbia). This development and 
flight success of SEADS clearly demonstrates the applicability of the FADS concept to blunt entry vehicles 
across the speed range. Through proper design, implementation and calibration an air data system based 
on the FADS concept can provide highly accurate air data to accomplish entry technology research and/or 
active flight control. 

Subsequent to the introduction of the FADS concept for the Shuttle Orbiter, the FADS concept has been 
recommended and/or considered and/or utilized to acquire air data on many flight systems including high 
performance aircraft 23 ’ 24 (B-2, F-14, F-18, F-22), guided munitions, 25 and entry/hypersonic systems with 
various degrees of bluntness. These systems include the blunt Aeroassist Flight Experiment (AFE) 26,27 and 
Pathfinder (which was not manifested) and higher fineness ratio vehicles including HYFLITE, HYFLEX, 
X-31, X-33, X-34, X-37, and X-43. FADS was recently utilized on the Orion Pad Abort One flight test, and 
will also be implemented on the Exploration Flight Test One. References 28-37 provide additional details 
as to the implementation of FADS concepts for these various applications. The flights of the SEADS and 
the design and development work accomplished in support of the AFE FADS have a heritage that closely 
relates to the MSL configuration. This background, coupled with the demonstrated flexibility of the FADS 
concept adequately demonstrates the desirability and the applicability of FADS to MSL. 

The MEDLI/MEADS project installed seven pressure ports through the MSL PICA heat shield at strate- 
gic locations to acquire heat shield surface pressure data during the atmospheric entry phase at Mars. The 
MEADS pressure ports are located on the MSL heat shield as shown in Figure 1(a). The PICA tile layout is 
also shown in this figure along with the predicted flow streamlines over the surface of the Thermal Protection 
System (TPS). All of the pressure ports are located a minimum of 3-inches from the PICA tile seams. At 
each of the pressure port locations there is an independent pressure measurement system installed on the 
internal surface of the heat shield support structure. 
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(a) MEDLI/MEADS Pressure Port Locations 


(b) MEDLI/MEADS Tube Configuration 


Figure 1: MEDLI/MEADS Geometry 


Pressure ports PI and P2 are located in the stagnation region to provide a nearly direct measurement of 
the total pressure at the high Mach regime. Ports P3, P4, and P5 lie on the spherical cap and are placed in 
order to take advantage of the simple geometry for angle of attack measurements. Additionally, P4, located 
at the geometric center, provides a nearly direct total pressure measurement at the low Mach regime prior to 
parachute deployment. The final two ports are located in the horizontal plane of symmetry, approximately 
1.0 meters off of the centerline. The ports P6 and P7 provide the off-axis measurements needed to estimate 
the angle of sideslip. The pressure ports are connected to pressure transducers via a tube system illustrated 
in Figure 1(b). The measured pressures are sampled at a rate of 8 Hz during entry. 

This paper provides an overview of the data processing algorithms utilized for MSL/MEADS trajectory 
reconstruction and the results of the reconstruction. These techniques consist of three semi-independent 
methods that utilize different portions of the entire data set. An overview of the raw data is provided, along 
with an assessment of data quality and comparison with simulation results. Techniques used to preprocess 
the data are described, which includes editing out suspect data and smoothing to reduce noise. The MISP 
sensors and associated modeling and reconstruction techniques are documented in Refs. 38-40 and are not 
discussed further in this paper. 

The remainder of the paper is organized as follows. Section II provides an overivew of the EDL timeline. 
Section III discusses the trajectory estimation techniques that have been implemented for MSL trajectory 
reconstruction. EDL data processing and trajectory reconstruction results are described in Section IV, and 
conclusions are provided in Section V. 

II. Mars Science Laboratory Entry, Descent, and Landing Overview 

Figure 2 shows a timeline of the different EDL events. 41 Note that the event times correspond to the 
nominal trajectory. Actual times differed from the nominal event times due to dispersions in the entry time 
and atmospheric uncertainties. EDL consists of six major segments: Exo-Atmospheric, Entry, Parachute 
Descent, Powered Descent, Sky Crane, and Fly Away. 42 The Exo-Atmospheric segment begins once the 
cruise stage separation command is sent. Once the cruise stage separates, Guidance, Navigation, and Control 
(GNC) is enabled. Once enabled, the entry body is despun and turned to its entry attitude. Then, the two 
75-kg Cruise Balance Masses (CBMs) are jettisoned to enable aerodynamic lift. 

The Entry segment starts with the vehicle at the Entry Interface Point (EIP) defined at 3522.2 km from 
the center of Mars, approximately 540 seconds after cruise stage separation. During the Entry segment, the 
vehicle goes through peak heating and peak deceleration, the Reaction Control System (RCS) controls the 
lift vector to achieve the desired down-range and cross-range target. Just prior to parachute deployment, 
six 25-kg Entry Balance Masses (EBMs) are jettisoned to eliminate lift and the vehicle rolls to point the 


3 of 25 


American Institute of Aeronautics and Astronautics 





Cruise Stage Separation 


Outside the Martian Atmosphere 




Cruise Balance Mass Jettison 
Turn to Entry Attitude 


a 


Entry Interface 

Peak Heating 

V 

V 

V 


Entry 


Peak Deceleration 

Heading Alignment 




Deploy Parachute 



(a) Cruise Stage Separation to Parachute Deployment 


(b) Parachute Deployment to Touchdown 


Figure 2: MSL Entry, Descent, and Landing 


Terminal Descent Sensor (TDS) to the ground. This maneuver is called the Straighten Up and Fly Right 
maneuver (SUFR) or “Victory” roll. Figure 3 shows the nominal trajectory over the time period from entry 
interface to parachute deployment. 

The Parachute descent segment starts with the parachute deployment triggered once the vehicle reached 
Mach 1.7. RCS wrist mode damping is active 10 s after parachute deployment to reduce wrist mode os- 
cillations. Once the vehicle achieves a speed of Mach 0.7, the heat shield is jettisoned and the TDS starts 
acquiring the ground. Note that the MEDLI instruments are powered off 10 seconds prior to heat shield 
jettison. The command to jettison the backsliell and the parachute is issued at an altitude of 1.6 km and 
at a velocity of approximately 79 m/s. Just before backshell separation, the Mars Landing Engines (MLEs) 
are primed in preparation for the start of the powered descent segment. 

The Powered Descent segment begins at backshell separation. During powered descent, eight indepen- 
dently throttleable MLEs are actuated, initially to execute a divert maneuver for backshell avoidance which 
brings the vehicle to vertical flight at a descent rate of 32 m/s. Once vertical flight is achieved, a descent 
at constant velocity to adjust for altitude error at backshell separation starts. This constant velocity phase 
is followed by a constant deceleration phase, which reducs the vehicle’s speed to 0.75 m/s in preparation 
for the sky crane segment. At this time, the four inboard MLEs are throttled down to near shutdown (1%) 
while the four remaining MLEs were throttled at 50%. 

The Sky Crane segment starts following issuance of the rover separation command, which occurrs at an 
altitude of approximately 18.6 m. The rover is lowered to 7.5 m below the descent stage. Then, the descent 
stage continues to descend until post-touchdown is detected. 

The Fly Away segment starts after touchdown is sensed. Once the descent stage stops it’s vertical motion, 
the bridle and electrical umbilical devices are cut and two of the MLE engines are throttled up to 100% while 
the other two engines are at slightly less than 100%. This causes the descent stage to pitch to 45 degrees. 
Once the turn maneuver is completed, all four engines are throttled up to 100%. Constant thrust is applied 
to ensure the descent stage impacts the surface at least 300 m from the landing point. 

III. Trajectory Reconstruction Methods 

This section provides a detailed description of various trajectory reconstruction methods that have been 
implemented for MSL EDL reconstruction. These reconstruction methods include three semi-independent 
methods that each make use of a different set of data for the reconstruction. The three semi-independent 
methods include an inertial reconstruction, that makes use of IMU measured accelerations and rates alone; an 
approach based on the vehicle aerodynamic database and sensed acceleration measurements; and a method 
based on the surface pressure measurements from the MEADS transducers. The first two of these methods 
have historically been the primary approaches to EDL reconstructions, for examples see References 43-48. 
The third method based on the surface pressure measurements has heritage from the SEADS program, but is 
new for this application to Mars EDL reconstruction. One key difference between the SEADS and MEADS 
pressure algorithms is that the SEADS pressure models were based on modified Newtonian theory, calibrated 
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(d) Aerodynamic and Trajectory Angles 


Figure 3: MSL Nominal Trajectory: Entry Interface to Heatshield Separation 


to wind tunnel experiments, whereas the MEADS pressure models are based entirely on CFD solutions. The 
following subsections provide details on each of the reconstruction approaches. 


A. Inertial Reconstruction 


The IMU-based reconstruction algorithm implemented for MSL EDL reconstruction is based on direct nu- 
merical integration of the equations of motion, using initial conditions provided by orbit determination 
solutions and star tracker attitude updates prior to cruise stage separation. The equations of motion used 
in the reconstruction are based on Reference 49-53, listed below. 
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where r is the radius of the vehicle from the center of the planet, 9 is the longitude, (f> is the declination, 
u, v, w are the inertial velocity components in a topocentric frame (defined by the z-axis direction along 
the radius vector, toward the center of the planet, and the y-axis to the east) are the Euler parameters 
describing the attitude of the vehicle with respect to the topodetic frame (north-east-down). The quantity 
Q is the planetary rotation rate. 

The body axis sensed accelerations are transformed from the IMU location according to the relations 
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where a x , a y , and a z are the sensed accelerations at the IMU location, x mi y rn and z m are the position of 
the IMU with respect to the center of mass in the body frame, and u> x , u> y and ui z are the body axis sensed 
angular velocity components. 

These equations of motion are integrated using a fixed step 4th-order Runge-Kutta algorithm with inte- 
gration rate of 200 Hz to match the rate of the MSL IMU data. Once integrated, the IMU reconstruction 
algorithm provides estimates of the planet relative trajectory based on the initial conditions and sensed 
accelerations and angular rates. Other quantities of interest, such as angle of attack, are computed from 
the integrated state variables during an output processing step. For IMU-only reconstruction, these outputs 
assume zero winds. An assumed atmosphere profile can be superimposed on the inertial trajectory in order 
to compute estimates of atmosphere-relative quantities such as dynamic pressure and Mach number. 

B. Aerodatabase Reconstruction 

This section provides a summary of the approach used to generate estimates of atmospheric relative param- 
eters based on the vehicle aerodynamic database, combined with accelerometer measurements. The method 
implemented for MSL closely follows the approach taken with Mars Pathfinder, described in Reference 44. 
The method first reconstructs the atmospheric density and pressure profile, and then solves for the angles 
of attack and sideslip from ratios of sensed accelerations. The entire process is iterated until the estimates 
converge. 



Figure 4: Aerodatabase Reconstruction Flow Diagram 

The block diagram shown in Figure 4 presents a high level description of the aerodatabase reconstruction 
algorithm. The algorithm requires the use of two loops designed to converge on accurate estimates of angle 
of attack, angle of sideslip, Mach number, density, pressure, and temperature. The outer loop performs the 
calculations to generate the air data parameters using outputs from the MSL aerodynamic database. The 
inner loop exists inside of the Newton solver to update angle of attack and sideslip. Both loops rely on 
the aerodatabase, specifically for force coefficient outputs. The MSL aerodynamic database and associated 
uncertainties are described in Refs. 54 and 55. 

The algorithm begins with an initial estimate of angle of attack, angle of sideslip, and Mach number. For 
the purposes of this application, an inertial reconstruction of these parameters is computed using the inertial 
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reconstruction method described in the previous section. During the first pass through the outer loop of the 
reconstruction algorithm, the MSL aerodatabase is called using the inertial inputs to generate an initial set 
of aerodynamic force and moment coefficients. The axial force coefficient is used to compute density. With 
this density estimate, pressure, speed of sound, Mach number, and temperature estimates can be computed. 



Figure 5: Aerodatabase Coordinate System 


These updated air data states are used as inputs to the aerodatabase to produce an updated set of force 
and moment coefficients to begin the inner loop of the algorithm. The inner loop of the algorithm uses the 
Newton method to solve a system of non-linear equations, which are functions of angle of attack and sideslip, 
and are composed of the acceleration and force coefficient ratios. Newton’s equation provides an update to 
angle of attack and sideslip using a numerically computed Jacobian matrix. The Newton solver is applied in 
a convergence loop that ends when the update to angle of attack and sideslip reaches a value smaller than the 
pre-defined threshold. (For this analysis, the convergence threshold was set to 10 -6 ). Once the loop ends, 
the updated angle of attack and sideslip are returned to the outer loop and used to query the aerodatabase 
in order to update the air data estimates. The outer loop and inner loop continue to run until both loops 
have converged. At this point, final estimates of angle of attack, angle of sideslip, pressure, density, speed 
of sound, Mach number and temperature have been reached. Mathematical details of the algorithm are as 
follows. 

First, the density is computed from the measured axial force and predicted axial force coefficient from 
the aerodatabase by means of 


2ma x 
P= SC A V 2 


( 7 ) 


where m is the vehicle mass, a x is the axial acceleration of the vehicle at the center of mass, V is the 
atmospheric relative velocity, S is the vehicle reference area, and Ca is the axial force coefficient from the 
aerodatabase. Once the density is obtained, the pressure can be computed from numerical integration of the 
hydrostatic equation, which is given by dp/dho = — pg , where p is the static pressure, g is the gravitational 
acceleration, and he is the radial altitude. This equation can be integrated along the trajectory using an 
Euler discretization of the form 

Pk = Pk-l - Pk9k ( h Gk - /iGfc-i) (8) 
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where k represents the time index. A simple gravity model can be used in the integration to compute 
acceleration as a function of altitude, given by 


9 = 9 o 


TD \ ^ 

R + ha) 


(9) 


where R is the planetary radius and go is the reference gravitational acceleration at Hq = 0. Once the density 
and pressure are obtained, the Mach number can be computed using the relationship M = V/ \fjpjp- 

Next, the angles of attack and sideslip are computed based on ratios of accelerations combined with a 
Newton solver. The approach utilizes a comparison between the ratio of accelerations and the corresponding 
ratio of force coefficients. Ratios used for the estimation of the angles of attack and sideslip are given by 
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where a x , a y , and a z are the body axis accelerations at the center of mass, and Ca, Cy, and Cn are the 
axial, side, and normal force coefficients from the aerodatabase. These ratios can be rearranged to form 
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(12) 


Next, a root finding algorithm is applied to solve for the angle of attack and sideslip that satisfy Equa- 
tion 12. For MSL reconstruction, a Newton method has been applied to solve Equation 12. The algorithm 
begins with an initial estimate of a and (3 based on the inertial reconstruction and then iterates until 
convergence using 

T+i } = { 7- } - 3 (a " A)_1 G (ai! fji} (13) 

where i denotes the iteration counter and J is the Jacobian matrix given by 


dgi/da dgi/d/3 
dg 2 /da dg 2 /dj3 


(14) 


which is evaluated numerically using perturbation values of 1 degree. 

At each time point, the Newton solver must iterate until the roots converge upon the final estimate of 
angle of attack and sideslip. Convergence is determined to have occurred when the difference between two 
consecutive iterations of the algorithm has fallen below a desired tolerance or when a maximum number of 
iterations has occurred. 

A key advantage of this reconstruction approach is that the estimates of the angles of attack and sideslip 
are wind-relative since the estimates are based on the aerodatabase. It is conceivable that the two different 
measures of the aerodynamic flow angles could be used to identify winds, in a manner similar that proposed 
in Reference 59; however, that approach suffers from a serious pitfall for the application to MSL so is not 
implemented here. The next section will elaborate further on issue of wind estimation. 


C. Pressure-Based Reconstruction 

The basic MEADS atmospheric state estimation procedure is essentially similar to the SEADS state esti- 
mation described in Reference 13, a point-wise least squares fit of the surface pressure measurements to the 
modeled pressure distribution to determine minimum-variance estimates of the atmospheric states. When 
combined with an Inertial Measurement Unit (IMU), the MEADS measurements can also be used to generate 
wind and density estimates. This section describes the techniques used to estimate these quantities, first 
providing details on the fore body surface pressure models and subsequently the basic MEADS estimation 
technique to compute the atmospheric state variables from the surface pressure measurements. 
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Figure 6: MEADS CFD Grid and Port Pressure 


1. Pressure Modeling 

A CFD-based table lookup model was developed for analysis of the MSL air data system. This CFD 
database is based on supersonic and hypersonic regime pressure distributions from thin layer Navier-Stokes 
solutions generated using the Langley Aerothermal Upwind Relaxation Algorithm (LAURA). The relaxation 
of vibrational temperature of the CO 2 molecule is simulated using the Camac 56 model. Forebody solutions 
were obtained using a 7 block singularity-free grid. 

The supplied CFD database was converted from absolute pressure to pressure coefficient in order to 
allow the pressure distribution to scale with different trajectories and to enable estimation of freestream 
atmospheric properties such as the static pressure. The data was re-interpolated from the 7 block grid to 
a single-zone grid of clock and cone angles. This alteration simplifies the interpolation of surface conditions 
without complications of search routines required to handle multiple zones. 

The CFD grid points are shown as a function of total angle of attack and Mach number in Figure 6(a) 
along with the nominal reference trajectory. The pressure profiles for each MEADS port are shown in 
Figure 6(b). Note that each CFD grid point consists of a full surface pressure distribution solution, with 
37 clock angles in uniform 5 deg increments and 61 cone angles with non-uniform increments. These grid 
points can be interpolated as needed to provide estimates of the pressure distribution at any point on the 
aeroshell. 


2. Atmospheric State Estimation 

By defining the atmospheric state vector as x , where x = [cr, /?, q,p s ] T , P is the vector of N observed surface 
pressures, and the CFD-based pressure model as h(x,t ), then the pressure measurement model can be 
written compactly as 

p = h(x,t k ) + e (15) 

where e is the vector of pressure measurement errors and t k is the measurement sampling time. The pressure 
measurement model can be approximated by means of the truncated series expansion 


p « h (x, tk) + H(tk) {x - x) + e 


(16) 


where x is some reference state and H(tk ) is the Jacobian matrix 


H(t k ) 


dh' 


dx 


x=x 


(17) 


which is also sometimes known as the measurement sensitivity matrix. The components of this matrix are 
shown in Figure 7 along the nominal trajectory. These sensitivities illustrate which ports provide the most 
information on a particular parameter, for instance the most leeward port 5 provides the most information 
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Figure 7: Pressure Sensitivities 


on the angle of attack, ports 6 and 7 provide the most information on sideslip, and ports 1 and 2 provide 
the most information on dynamic pressure. 

The aerodynamic state estimation problem can be reduced to a linear regression problem of the form 

y = Hx + e (18) 

where y = p — h (x, tp) + Hx. By virtue of the Gauss-Markov theorem, the best linear unbiased estimate 
of x is the weighted least-squares solution, 58 

x= (H T R- 1 H') 1 H T R 1 y (19) 

where R is the pressure measurement error covariance matrix. Since the original relationship between the 
states and the measurement in Eq. (15) is nonlinear, the estimation scheme can be iterated until convergence 
by successively replacing x by x , using the converged state estimate from the previous sampling time (tk- i) 
as the initial reference state x. The state estimate error covariance matrix P can then be computed from 

P = (iT 7 # -1 #) 1 (20) 

The algorithm can also be iterated globally in order to estimate and remove systematic errors that result 
from the combination of transducer and CFD distribution uncertainties. After each pass through the pressure 
data, the residuals are computed as r = p — p where p are the measured pressures and p are the modeled 
pressures based on the converged state. The time history of the residuals for each port are then fitted with 
a second-order polynomial of the form 

n= cii + biPi + Cip \ ( 21 ) 
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where r* is the residual for the ?’th pressure port. The coefficients of the polynomial are determined using 
a least-squares fit to the residual time history, and are then used to correct the measurements for the next 
global iteration. 

The least-squares solution described in this section has several important properties. First, all pressure 
measurements are processed simultaneously to determine the atmospheric state variable estimate, opposed to 
the so-called “triples” algorithms, which use only a subset of available ports to determine the state variable 
estimates individually in sequence. Using all measurements simultaneously is advantageous in the sense that 
every port contributes to the state estimate rather than a restricted set, which generally leads to better 
performance, 57 and is also fault tolerant in the sense that a failed port can simply be removed from the 
solution that would otherwise ruin a triples algorithm. 

The aerodynamic state vector formulation is not unique. In particular, any two variables from the set of 
p s , pt , q , R, and M can be used in the estimator in addition to a and /?. The remaining state variables not 
used directly in the estimator can be calculated as parameters, having first estimated the aerodynamic state 
variables. Note that the SEADS algorithm 13 made use of p s and pt in the estimator formulation. 

The quantity R listed above is the pressure ratio. The pressure ratio is a function of the Mach number, 
which takes the form 13 
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(22) 


The pressure ratio is shown as a function of Mach number for 7 = 1.335 in Fig. 8(a). 

The pressure measurement weighting as needed for the MEADS-based pressure reconstruction consists 
of system error models of transducer characterizations from thermal vacuum chamber calibrations, thermo- 
couple measurement errors such that an inaccurate temperature is used in the calibration database, system 
noise and quantization, time tag errors and sampling delays, pressure path leaks, pneumatic lag, and ther- 
mal transpiration. Additionally, the MEADS weighting includes a model of the CFD uncertainties, which 
consists of deformation, ablation, outer mold line changes, grid refinement, protuberances, and port location 
uncertainty. These error models are summarized in Refs. 60-62. 


3. Atmosphere Reconstruction 


Other atmospheric relative parameters can be estimated from the MEADS pressure solution in conjunction 
with the IMU measurements and velocity estimates, including the atmospheric density and the basic vehicle 
force and moment aerodynamic coefficients. First, The freestream atmospheric density follows from the 
determination of the dynamic pressure from 



(23) 


where V is the IMU-based planet-relative velocity magnitude. The uncertainty in the density estimate can 
be computed as 
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where a q is the dynamic pressure uncertainty from the least-squares state covariance matrix, and ay is the 
relative velocity uncertainty computed from linear covariance analysis of the IMU-based reconstruction. 

Once computed, the density profile can be integrated using the hydrostatic equation method shown in 
Eq. (8) to produce another estimate of the freestream static pressure. Then, the density and static pressure 
estimates can be used to compute an estimate of atmospheric temperature from the equation of state. 

It is possible to estimate winds by having both a wind-relative estimate of the aerodynamic flow angles 
from the pressure-based solution and a planet-relative estimate from the inertial reconstruction, using an 
approach based on Reference 59. In this approach, the wind components are determined by solving a 
nonlinear system of algebraic equations that arise by equating the winds to a function of the difference 
between the wind-relative and planet-relative aerodynamic angles. This method was implemented and tested 
for application to MSL in Reference 60, but a drawback of the method is that it is singular for zero flight path 
angle, which presents a serious limitation for application to MSL. Therefore, the method is not proposed for 
use in MSL EDL data reduction. 
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4- Aerodynamic Reconstruction 


Aerodynamic force and moment coefficients can be estimated by combining the atmospheric state estimates 
with the sensed accelerations and angular rates from the IMU. The force coefficients can be computed from 
the equations 


Ca 
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ma y 

Sq 

ma z 
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(26) 
(27) 


where C A , Cy, and Cn are the axial, side, and normal force coefficients, respectively. Here, a x , a y , and 
a z are the body axis sensed accelerations at the vehicle center of mass, m is the vehicle mass, and S is the 
aerodynamic reference area. 

Since the mass, acceleration, and dynamic pressure uncertainties are uncorrelated, simple expressions can 
be obtained for the uncertainty in the derived axial, normal, and side aerodynamic force coefficients. The 
expression for the axial force coefficient uncertainty is given in the following equation. 
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Similar expressions can be obtained for the side and normal force coefficients by substituting the appro- 
priate acceleration axis subscript. 

A similar approach can be used to determine moment coefficients by making use of angular velocity 
measurements from the IMU and the estimated spacecraft inertia. 


5. Inertial Measurement Unit Aided Estimation Algorithm 

As noted in Ref. 13-14, for small pressure ratios (or, equivalently, high Mach) the product RM 2 is nearly 
constant. This trend is illustrated in Fig. 8(b) for 7 = 1.335. The implications of this trend are twofold. 
First, any slight change in static pressure has a large impact on the Mach number, which implies that the 
Mach number is difficult to estimate for low pressure ratios. Secondly, it implies that the dynamic pressure 
can be computed accurately directly from the total pressure by substituting the limiting value of the quantity 
RM 2 . Therefore, errors in static pressure or Mach number do not contribute to errors in dynamic pressure 
estimates in the high hypersonic regime. 

Due to the asymptotic behavior of the pressure ratio as a function of Mach number, some numerical 
difficulties in the FADS estimation algorithm can arise for high Mach number conditions. For high Mach 
cases it can be advantageous to consider a special case of estimators that use a fixed pressure ratio and form 
state estimates using only three aerodynamic variables, e.g. pt, a, and /3, with dynamic pressure computed 
directly from p t assuming the limiting value of the quantity RM 2 . 63 These estimators are generally more 
consistent but become biased as the true pressure ratio drifts away from the assumed value. If the vehicle 
is instrumented with an IMU then it is also possible to schedule the assumed pressure ratio as a function of 
the inertial Mach number to reduce the growth of the bias over time. Alternately, the bias can be calibrated 
and subtracted out using preflight simulation cases. 

For application to MSL, a new approach is taken in which the IMU-derived Mach number is used as 
a pseudo-observation in the least-squares algorithm. In this approach, another equation is added to the 
observation vector, leading to a system of equations given by 

2 = h (x, tk) + e (30) 

where 
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Figure 8 : Pressure Ratio vs. Mach for 7 = 1.335 
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Figure 9: Mach Number Measurement Sensitivities 
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Following Sec. III.C.2, the weighted least-squares solution is given by 

/ ~ t — 1~\ -1 ~ t — 1 

x = (H R HJ H R y 


(32) 


(33) 


where y = z — h (x, t k ) + Hx , H — [dh/dx\ x =xi and R is the measurement covariance matrix that is aug- 
mented with an uncertainty estimate for the Mach number pseudo-observation. The additional components 
of the measurement sensitivity matrix for the IMU-aiding algorithm are shown in Fig. 9. 

Eq. (33) is iterated at each observation time until the state estimate converges (or a maximum number 
of iterations has been reached). The covariance matrix for the state estimate is given by 

P = [h T R 1 H^j 1 (34) 

To apply the IMU aiding to the least-squares FADS algorithm, an estimate of the speed of sound is needed 
to compute Mach number from the IMU-derived velocity. This speed of sound estimate can be based on the 
best available prior atmosphere knowledge. Since the FADS algorithm computes a speed of sound estimate 
based on the computed density and static pressure, the algorithm can be iterated globally by updating the 
speed of sound estimate on each pass, based on the derived speed of sound from the previous iteration. 
Systematic errors are also estimated on each global iteration of the algorithm as described in Sec. III.C.2. 
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IV. Mars Science Laboratory Entry, Descent, and Landing Data Analysis 

A. Data Quality and Preprocessing 

The following sections provide details of the various EDL sensor data that were acquired for trajectory and 
atmosphere reconstruction. These data consist of IMU and pressure sensor data. Also, atmospheric models 
required for the pure inertial trajectory reconstruction are described. 

1. Inertial Measurement Unit Data 

Data from the onboard IMU, consisting of AVs and Ads were acquired at a rate of 200 Hz. These velocity 
and attitude increments were converted into acceleration and angular rate by dividing by the time step for 
each measurement. A Tustin bilinear filter was used to smooth the data for use in the aerodatabase and 
MEADS aerodynamic reconstructions. Unfiltered data was integrated for the pure-IMU reconstruction. The 
filtered and unfiltered IMU data are shown in Fig. 10 for the period from entry interface (540 seconds) 
to touchdown. Note that considerable noise due to structural vibrations occur during the bank reversal 
maneuvers where the RCS jets are active. The initial conditions for the integration are based on the orbit 
determination (OD) number 229, which provides the position and velocity of the vehicle in the Earth Mean 
Equator of J2000 (EMEJ2000) coordinate system at a time of approximately 550 seconds prior to entry 
interface. 66,6 ' The attitude initialization was based on the on-board navigation filter that performed a star 
tracker alignment prior to cruise stage separation. This attitude solution was then propagated to the same 
initial time as the OD229 solution by numerical integration. 

2. Pressure Data 

The MEADS pressure measurement data were acquired at an 8 Hz sample rate during EDL. For EDL 
reconstruction, the data were first converted into engineering units, after which an in-flight zero was applied 
to correct for transducer thermal drift. Subsequently, outlying data points due to the entry ballast mass 
ejections were edited out and filled in with a first-order polynomial fit, and a 1 Hz optimal Fourier filter 
was applied to smooth the data. The pressures were then interpolated to the Port 4 time tag in order to 
generate pressure samples at a common 8 Hz rate. The resulting pressures are compared to pre-flight nominal 
trajectory predictions based on OD 213 in Figs. 11- 13. Fig. 11 shows comparisons for the stagnation region, 
Ports 1 and 2. Fig. 12 shows the comparisons for the ports along the spherical cap, Ports 3-5, and Fig. 13 
shows comparisons for Ports 6 and 7, which provide the primary measurement of the sideslip angle. These 
figures illustrate that the sensor readings during EDL are reasonable, and that all transducers are functional. 
Data quality during EDL was good for the majority of the trajectory. Some noise was introduced during 
pyro events associated with mass ejections during the reorientation to zero angle of attack in preparation 
for parachute deployment. Vibrations induced by the parachute mortar fire event were large, which further 
reduced the MEADS data quality. This vibrational noise, coupled with the rapid decrease in dynamic 
pressure following parachute deployment, essentially limited the last valid MEADS measurement to the 
instant prior to parachute mortar fire. Other aspects of MEDLI hardware performance during EDL are 
described in Ref. 64. 

3. Atmospheric Models 

Atmospheric models along the MSL EDL trajectory were generated from preflight mesoscale models, de- 
scribed in Ref. 65, tuned to match surface pressure measurements of 695 Pa from Curiosity, which were 
obtained after landing. The modeled atmosphere is an average of two such mesoscale models, namely 
the Mars Mesoscale Model 5 (MMM5) and the Mars Regional Atmospheric Modeling System (MRAMS). 
These models provide data up to 50 km altitude, which can be extrapolated to higher altitudes as needed. 
This combined model was queried along the IMU-navigated trajectory at time increments of 0.5 seconds 
to generate pressure, density, and temperature profiles. These profiles were then interpolated by altitude 
in order to provide atmospheric-relative estimates (Mach number and dynamic pressure) for the pure IMU 
reconstruction. 
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Figure 10: Sensed Accelerations and Angular Rates 


B. Reconstruction Results 

Results of the three independent reconstructions are described in the following subsections. These sections 
describe the reconstructions of the aerodynamic state, free stream atmosphere, and vehicle aerodynamics. 

1. Aerodynamic State Reconstruction 

The aerodynamic state reconstructions are shown in Fig. 14. Fig. 14(a) shows results for the angle of attack, 
Fig. 14(b) shows results for the angle of sideslip, Fig. 14(c) shows results for dynamic pressure, and Fig. 14(d) 
shows results for the Mach number estimate. In each of these plots, the pure-IMU reconstruction is shown in 
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(a) Port 1 (b) Port 2 


Figure 11: MEADS Port 1 and 2 Pressures 




(a) Port 3 (b) Port 4 



(c) Port 5 

Figure 12: MEADS Port 3-5 Pressures 


blue, the aerodatabase reconstruction is shown in red, and the IMU-aided MEADS pressure reconstruction is 
shown in cyan. The three reconstructions give fairly consistent results. The hypersonic trim angle predicted 
by the pure-IMU reconstruction closely matches the predicted trim angle based on the aerodatabase, which 
is shown in the black curve in Fig. 14(a). There are differences in angle of attack on the order of 0.4 degrees 
between the inertial, MEADS, and ADB reconstructions in the hypersonic regime. The three reconstructions 
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(a) Port 6 (b) Port 7 


Figure 13: MEADS Port 6 and 7 Pressures 
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Figure 14: Aerodynamic State Reconstruction 


are in better agreement in the supersonic regime. 

2. Atmospheric Reconstruction 

The atmospheric reconstructions are shown in Fig. 15. Fig. 15(a-b) show results for the estimate of density, 
Fig. 15(c-d) show results of the static pressure estimate, and Fig. 15(e-f) show the results of the temperature 
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(a) Density vs. Time 


(b) Density vs. Altitude 



(c) Pressure vs. Time 



(d) Pressure vs. Altitude 



(e) Temperature vs. Time 
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Figure 15: Atmosphere Reconstruction 


reconstruction. In each case, the left plot shows the results as functions of time along the trajectory whereas 
the right plot shows the results as a function of altitude. Note that the IMU reconstruction here refers 
to the mesoscale model queried along the reconstructed altitude profile. Two MEADS reconstructions of 
static pressure are shown in Fig. 15(c-d). Here, the cyan curve shows the static pressure estimate that is 
produced directly from the least-squares processing algorithm, whereas the dashed green curve shows the 
static pressure estimate that is produced from a hydrostatic integration of the derived density profile. These 
two reconstructions are in close agreement, with maximum differences on the order of 0.02%. The initial 
condition for the hydrostatic integration was obtained from the combined mesoscale atmosphere model at 
an altitude of 66 km. The three reconstructions are in fairly good agreement until an altitude of 13-14 km. 
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Recall that the MEADS derived density assumes that the atmospheric winds are zero relative to the 
rotation of the planet. As the vehicle decelerates, any non-zero wind becomes larger in proportion to the 
IMU-based planet relative velocity, and these winds produce an error in the density estimate. A question 
that naturally arises is when and at what altitude do these potential errors become important to the recon- 
struction. To answer this question, a simple wind sensitivity study was conducted in which a 20 m/s tail or 
head wind (represented by a pure east/west wind component) was included in the estimated relative velocity 
based on the IMU integration. The results are shown in Fig. 16(a) and (b), where the results are shown as 
functions of time and altitude, respectively. The results of the nominal zero wind case are also shown for 
comparison. The differences between the reconstructions starts to become apparent at a time of roughly 
640-650 seconds, corresponding to an altitude of 13-14 km. Interestingly, this is the same altitude range 
where the nominal MEADS reconstruction begins to differ from the mesoscale model. These results suggests 
that the differences between the MEADS reconstruction and the mesoscale model could be reconciled with 
a reasonable tail wind. 




(a) MEADS Wind Sensitivity vs. Time 


(b) MEADS Wind Sensitivity vs. Altitude 


Figure 16: MEADS Wind Sensitivity 


3. Aerodynamic Reconstruction 
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Figure 17: Aerodynamic Reconstruction 


The MEADS-based reconstruction of the vehicle aerodynamics is shown in this section. Fig. 17(a) shows 
the axial force reconstruction and Fig. 17(b) shows the normal force reconstruction. In each case, the 
MEADS-based reconstruction is shown in the cyan curve, and the nominal aerodatabase queried at the 
MEADS-derived atmospheric flight condition is shown in the blue curve. The MEADS results indicate 
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a slightly higher drag than the nominal prediction. Complete analysis and model reconciliation of the 
aerodynamic reconstruction is beyond the scope of this paper. More detailed analysis of the MSL entry 
vehicle aerodynamic performance are provided in Ref. 68. 

4- MEADS Pressure Residuals 
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Figure 18: Port 1 and 2 Residuals 


This section shows the pressure residuals from the IMU-aided MEADS reconstruction. Recall that the 
first pass of the algorithm assumes an isothermal atmosphere with no systematic errors, and subsequent 
passes iterate to solve for the atmospheric conditions and to determine systematic errors. The pressure 
residuals are expected to reduce as a consequence of these iterations. The results based on the MEADS 
flight data are shown in Figs. 18-20. These plots show the residuals for Ports 1-2, 3-5, and 6-7, respectively. 
In each case, the residual based on the first iteration is shown in the blue curve and the residual based on the 
final converged solution is shown in cyan. Also, the 3cr uncertainty bands that include measurement system 
uncertainties as well as CFD pressure model uncertainties are shown. The left column of figures shows the 
residual in units of pressure whereas the right column shows the residuals as percent of reading. These 
plots indicate that the IMU-aided MEADS algorithm is successfully able to reduce the residuals and remove 
systematic error by iterating on the atmosphere and systematic error parameters. The resulting residuals are 
well within the expected 3cr uncertainty bands and roughly zero mean, which is another indication of good 
transducer performance during EDL. The Root Mean Square (RMS) residuals and estimates of systematic 
error parameters over the period of good quality data (roughly 590-775 seconds) are summarized in Table 1. 
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Figure 19: Port 3-5 Residuals 


21 of 25 


American Institute of Aeronautics and Astronautics 



(a) Port 6: Pascals 
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Figure 20: Port 6 and 7 Residuals 


Table 1: RMS Pressure Residuals and Systematic Error Estimates 


Port 

Residual RMS Error 


Systematic Error Estimates 


First Iteration (Pa) 

Final Iteration (Pa) 

Bias (Pa) 

Scale Factor 

Nonlinearity (1/Pa) 

1 

62.76 

5.14 

1.75 

-6.28E-03 

3.72E-08 

2 

91.31 

11.30 

-1.67 

-6.80E-03 

-3.23E-08 

3 

234.55 

14.35 

-1.77 

2.23E-02 

-9.10E-08 

4 

115.16 

18.70 

-25.51 

1.57E-02 

-1.59E-07 

5 

110.73 

19.41 

11.67 

-2.41E-02 

5.65E-07 

6 

80.68 

9.50 

0.56 

-6.05E-03 

-1.41E-07 

7 

52.09 

6.07 

-0.99 

-3.50E-03 

-1.08E-07 
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V. Conclusions 


This paper describes three complimentary reconstruction methods that have been implemented for Mars 
Science Laboratory entry, descent, and landing trajectory reconstruction. Each of these methods makes use 
of a different subset of the total available data to reconstruct portions of the trajectory. Specifically, the 
methods include a pure inertial reconstruction, a reconstruction based on the nominal aerodatabase, and a 
reconstruction based on measurements of the forebody surface pressure distribution. The first two of these 
methods are considered standard approaches for Mars entry, descent, and landing reconstruction, whereas 
the latter of these techniques is new to this application. The pressure-based reconstruction is enabled by 
instrumentation consisting of transducer measurements from seven pressure ports on the vehicle heatsliield. 
Such heatshielcl pressure data has not successfully been collected during Mars entry prior to Mars Science 
Laboratory. A novel inertial measurement unit-aided estimation approach is introduced in this paper for 
processing these pressure measurements, wherein the planet-relative velocity is utilized to help anchor the 
pressure-derived Mach number estimates, which otherwise are nearly unobservable with pressures alone. An 
important feature of the new pressure-based reconstruction technique is that it allows for an atmosphere 
reconstruction to be conducted without assumptions on the vehicle aerodynamics. This in turn allows for 
an aerodynamic reconstruction to be performed, which has not been possible for Mars entry reconstructions 
of the past. 

These algorithms are utilized to process flight data from Mars Science Laboratory’s atmospheric entry at 
Mars, which occurred on August 5th 2012. Assessments of the pressure measurements indicate that the data 
quality is good. Vibrations from pyrotechnic events effectively limit the last valid pressure measurement to 
the instant prior to parachute mortar fire. The three reconstruction methods are generally in good agreement. 
The pure inertial reconstruction of angle of attack is in close agreement to the preflight predicted trim angle. 
The reconstructed atmosphere profile from the pressure reconstruction is in reasonable agreement with 
atmospheric models that were tuned to match surface pressure measurements from the Curiosity rover until 
an altitude of approximately 13-14 km. In this altitude range the pressure-based reconstruction becomes 
sensitive to atmospheric winds, which could account for the drift away from the modeled atmosphere profile. 
The overall good data quality and general agreement between the reconstructions supports a future combined 
reconstruction involving Kalman filtering techniques to optimally blend all available data together into a 
single trajectory and atmosphere estimate. 
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